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Abstract 

A novel scenario for the emergence of geometry in random multitrace matrix models 
of a single hermitian matrix M with unitary U{N) invariance, i.e. without a kinetic term, 
is presented. In particular, the dimension of the emergent geometry is determined from 
the critical exponents of the disorder-to-uniform-ordered transition whereas the metric is 
determined from the Wigner semicircle law behavior of the eigenvalues distribution of the 
matrix M. If the uniform ordered phase is not sustained in the phase diagram then there 
is no emergent geometry in the multitrace matrix model. 

1 Introduction and Motivation 

The original motivation for this work is the theory of noncommutative <I>^ which we now 
briefly describe. A scalar phi-four theory on a non-degenerate noncommutative Euclidean 
spacetime is a a three-parameter matrix model of the generic form 

S = Tr^ (oM AM + + cM^). (1.1) 

The Laplacian A captures precisely the underlying geometry, i.e. the metric, of the 

noncommutative Euclidean spacetime in the sense of m- This theory can be regularized 
non-perturbatively using N x N matrices in an almost obvious way, i.e. the Hilbert space 
H can be taken to be finite dimensional of size N. This theory exhibits the following 
three known phases: 

• The usual 2nd order Ising phase transition between disordered < M >= 0 and 

uniform ordered < M >~ 1 w phases. This appears for small values of c. This is the 

only transition observed in commutative phi-four, and thus it can be accessed in a 
small noncommutativity parameter expansion. 
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• A matrix transition between disordered < M >= 0 and non-uniform ordered < 

M >~ 7 phases with 7 ^ = Itv- This transition coincides, for very large values of c, 
with the 3rd order transition of the real quartic matrix model, i.e. the model with 
a = 0, which occurs at h = —2\/Nc. In terms of 6 = and c = cN~‘^ this 

reads 

b=-2y/l. (1.2) 

This is therefore a transition from a one-cut (disc) phase to a two-cut (annulus) 
phase [281130] . See also [291I3T] . 

• A transition between uniform ordered < M >~ Itv and non-uniform ordered < 
M >~ 7 phases. The non-uniform phase, in which translational/rotational invari¬ 
ance is spontaneously broken, is absent in the commutative theory. The non-uniform 
phase is essentially the stripe phase observed originally on Moyal-Weyl spaces in m- 

Thus, the uniform ordered phase < $ >~ In is stable in the theory dni). This fact is in 
contrast with the case of the real quartic matrix model V = TvnibM"^ -|- cM^) in which 
this solution becomes unstable for all values of the couplings. The source of this stability 
is obviously the addition of the kinetic term to the action. 

The non-uniform ordered phase m is a full blown nonperturbative manifestation of 
the perturbative UV-IR mixing effect [26] which is due to the underlying highly non-local 
matrix degrees of freedom of the noncommutative scalar field. 

The above picture of the phase diagram holds for noncommutative phi-four in any 
dimension, and the three phases are all stable, and are expected to meet at a triple 
point. The phase structure in four dimensions was discussed using the Hartree-Fock 
approximation in [T] and studied by means of the Monte Carlo method, employing the 
fuzzy torus |37| as regulator, in [2]. 

In two dimensions the noncommutative phi-four theory is renormalizable [3|. The 
regularized theory on the fuzzy sphere W is given by the action (HD with a finite 
dimensional Hilbert space H of size N and a Laplacian A = [La, [La, ..]] where La are the 
generators of SU{2) in the IRR of spin {N — l)/2. 

The above phase structure was confirmed in two dimensions by means of Monte Carlo 
simulations on the fuzzy sphere in pnoj- Indeed, fuzzy scalar phi-four theory enjoys three 
stable phases: i) disordered (symmetric, one-cut, disk) phase, ii) uniform ordered (Ising, 
broken, asymmetric one-cut) phase and hi) non-uniform ordered (matrix, stripe, two-cut, 
annulus) phase. The phase diagram is shown on the two graphs of figure ([T]) which were 
generated using the Metropolis algorithm. 

The problem of the phase structure of fuzzy phi-four was also studied by means of 
the Monte Carlo method in [nHiai25|. The analytic derivation of the phase diagram of 
noncommutative phi-four on the fuzzy sphere was attempted in [T5l[20l[2n[32] . 

The related problem of Monte Carlo simulation of noncommutative phi-four on the 
fuzzy torus, and the fuzzy disc was considered in la, m, and [22| respectively. For a 
recent study see [23] . 

In [ 25 ] the phase diagram of fuzzy phi-four theory was computed by Monte Carlo 
sampling of the eigenvalues Aj of the scalar field M. This was possible by coupling the 
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scalar field M to a 17(1) gauge field Xa on the fuzzy sphere which then allowed us, 
by employing the U{N) gauge symmetry, to reduce scalar phi-four theory to only its 
eigenvalues. The pure gauge term is such that the gauge field Xa is fluctuating around 
Tfa = La. 

Another powerful method which allows us to reduce noncommutative scalar phi-four 
theory to only its eigenvalues, without the additional dynamical gauge field, is the mul¬ 
titrace approach. The multitrace approach was initiated in [151116] . See also [23] for a 
review and an extension of this method to the noncommutative Moyal-Weyl plane. For 
an earlier approach see [32| and for a similar more non-perturbative approach see UZHSoj. 
The multitrace expansion is the analogue of the Hopping parameter expansion on the 
lattice in the sense that we perform a small kinetic term expansion, i.e. expanding in the 
parameter a of m, while treating the potential exactly. This should be contrasted with 
the small interaction expansion of the usual perturbation theory. The effective action 
obtained in this approach is a matrix model which can be expressed solely in terms of the 
eigenvalues A* and which, on general grounds, can only be a function of the combinations 
T 2 n oc ~ To the lowest non-trivial order we get an effective action of the 

form millZllll] 


5efr = 5](6A2 + cAf)--^ln(Ai-A,)' 




+ 


aN 




— '<^ 2,1 - Aj)^ - - Aj)^] -h 








(1.3) 


The logarithmic potential arises from the Vandermonde determinant, i.e. from diagonal- 
ization. The coefficients U 2 ,i, and ^ 2,2 are given by U 2 ,i = +1 , ^ 4^1 = 0 , ^ 2,2 = 1/8. 

Furthermore, it is not difficult to convince ourselves that the above action is a multitrace 
matrix model since it can be expressed in terms of various moments = TrM"^ of the 
matrix M. 

The original multitrace matrix model written down m comes with different values of 
u’s and therefore, in the commutative limit N —> 00 , it corresponds to a phi-four theory 
on the sphere modulo multi-integral terms. 

Since these multitrace matrix models depend only on N independent eigenvalues their 
Monte Carlo sampling by means of the Metropolis algorithm does not suffer from any 
ergodic problem. The phase diagrams of these models obtained in Monte Carlo simulations 
will be reported elsewhere. 

The remainder of this article is organized as follows; 

1. Section 2: We describe our proposal for how fuzzy geometry can emerge in generic 
multitrace matrix models. 

2. Section 3: We apply our proposal to an explicit example. We will show that if the 
multitrace matrix model under consideration does not sustain the uniform ordered 
phase then there is no emergent geometry. On the other hand, if the uniform ordered 
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phase is sustained then there is an underlying or emergent geometry. In particular, 
we will show how 

• i) to determine the dimension from the critical exponents of the uniform-to- 
disordered (Ising) phase transition, and how 

• ii) to determine the metric (Laplacian, propagator) from the Wigner semicircle 
law behavior of the eigenvalues distribution of the matrix M. 

3. Section 4: We conclude by giving a straightforward generalization to fuzzy CP” and 
fuzzy T”. 


2 The Proposal 


We start with a general multitrace matrix model rewritten in terms of the moments 
TrM” with generic parameters B, C, D, B', C', D', ^ ,...as 

1 2 


V = BTrM'^ + CTrM^ + D 


TrM^ 


+ b'[T rMf+ C'TrMTrM^ + d'[T rM)^ + XTrM‘^{TrMf + .... (2.1) 


This action includes the noncommutative phi-four model on the fuzzy sphere (jl.3l) and 
the multitrace matrix model of |15j as special cases. It also includes as special cases the 
multitrace matrix models obtained by expanding the kinetic term on i) fuzzy CP” |16p36j . 
on ii) Moyal-Weyl spaces with and without the harmonic oscillator term [24], and on hi) 
fuzzy tori m- 

The phase diagram of the action (12.ip will generically contain the matrix one-cut-to- 
two-cut transition line separating the two stables phases of disorder and non-uniform- 
order. However, the uniform ordered phase will typically be unstable as in the case of the 
real quartic matrix model 

V = BTrM^ + CTrM^. (2.2) 


Our proposal goes as follows. We can check for a possible emergence of geometry in the 
multitrace matrix model (12.Ih by following the three steps: 

1. We compute the phase diagram of the model (j2.ip . If the uniform ordered phase 
remains unstable as in the case of the real quartic matrix model (12.2p then there is 
no geometry and the model is just a trivial deformation of (|2.2p . In the opposite case 
we claim that there is an underlying, i.e. emergent, geometry with a well defined 
dimension (step 2) and a well defined Laplacian/metric (step 3). This means that we 
can rewrite the multitrace matrix model, in the region of the phase diagram where 
the uniform ordered phase exists, in terms of a scalar function and a star product 
with a noncommutativity parameter 6 by finding the appropriate Weyl map. As 
a consequence, a small noncommutativity parameter expansion can be performed 
and the the limit 9 —> 0 can be taken. The disordered-to-uniform-ordered phase 
transition reduces therefore to the usual 2nd order Ising phase transition on the 
underlying geometry. 
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2. We compute the dimension of the underlying by computing the critical exponents 
of the disordered-to-uniform-ordered phase transition which, by universality, take 
specific values in each dimension. 

3. We compute the Laplacian by computing the free behavior of the propagator. This is 
done explicitly by computing the eigenvalues distribution of the matrix M in the free 
regime, small values of C, and comparing with the Wigner semicircle law behavior 
which must hold with a specific radius depending crucially on the kinetic term. 

3 Explicit Example: The Fuzzy Sphere 

3.1 Phase Diagram 

We consider as an example the multitrace matrix model of m which comes with the 
V values U 2 ,i = —1 j va,i = 3/2 , U 2,2 = 0. The action is given explicitly by 

V = BTrM^+ CTrM^ + D[TrM^]^ + B'(TrMf+ C'TrMTrM^. (3.1) 

The parameters B, B' and C' are constrained as D = 3A^/4, B' = '/N /2 and C' = —N. 
The phase diagram of this model is computed by means of Monte Carlo elsewhere. The 
result is shown on figure ([2]). The details of the corresponding non-trivial lengthy Monte 
Carlo calculation will be reported elsewhere. As desired we have three stables phases 
in this particular model meeting at a triple point. In other words, we have established 
that this multitrace matrix model sustains the uniform ordered phase which is the first 
requirement. 

3.2 Dimension from Critical Exponents 

The uniform ordered phase is also called the Ising phase precisely because we believe 
that the corresponding transition to the disordered phase is characterized by the universal 
critical exponents of the Ising model in two dimensions derived from the Onsager solution. 
These critical exponents are defined as usual by the following behavior 

m/N =< \TrM\ > /A^ ~ {B^ - Bf ~ 

C^/N'^ ^ {B- Bc)~^ ~ 

X =< iTrMp > - < \TrM\ {B - Bc)~^ ~ ~ 

i^\B- Bc\~'" ~ N. (3.2) 

There are in total six critical exponents, the above five plus the critical exponent 6 which 
controls the equation of state, but only two are truly independent because of the so-called 
scaling laws. The Onsager solution of the Ising model in two dimensions gives the following 
celebrated values |33] 


v = 1 , fi = 1/8 , 7 = 7/4 , a = 0 , 7 = 1/4 , <5 = 15. 
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(3.3) 


This fundamental result is very delicate to check explicitly in the Monte Carlo data. Since 
we must necessarily deal with the critical region we must face the two famous problems of 
finite size effects and critical slowing down. In this particular problem, the critical slowing 
down problem can be shown to start appearing in Monte Carlo simulations around iV > 60 
so we will keep below this value and employ very large statistics of the order of 2 ^*^ to 
avoid it. A more systematic solution to this problem is to employ the Wolf algorithm |34] 
which we do not attempt here. We simply employ here the ordinary Metropolis algorithm. 
The problem of finite size effects is also very serious for the measurement of the critical 
exponents since the above behavior is supposed to hold only for large N. This 

problem can be avoided by not including values of N less than 20 and thus below we will 
quote for completeness = 10 and = 15 data but, in most cases, we will not take 
them into account in the fitting. 

Since the Ising model appears from the theory for large values of the quartic cou¬ 
pling it is preferable to use values of C as large as possible. However, we are limited from 
above by the appearance of the different physics of the transition between the disordered 
and non-uniform-ordered phases around C = 1.5. Thus, we choose C = 1.0 which is 
relatively large but well established to be within the Ising transition with an extrapolated 
critical point around B = —3.07 (see below). The critical behavior of the magnetization, 
susceptibility and specific heat around the critical value of 7? = —3.10 is shown on fig¬ 
ure ([3|). We attach in table (HJ some data relevant for the computation of the critical 
exponents v, /3, 7 and a. The other critical exponents can be determined via scaling laws. 

The measurements of the critical exponents v, /3, 7 and a proceeds as follows: 

• Critical Point and The Critical Exponent ly: By plotting the critical point 
Be obtained for each N versus N (first and second columns of table ([T])) we get 
immediately both the N = 00 critical point and the critical exponent v. We obtain 
(see figure dH)) 

Be = -1.061(168).Ar-°'926(83) _ 3 074(6) ^ , 1 / = 0.926(83). (3.4) 

Also we obtain 

= -3.074(6). (3.5) 

This prediction for ly agrees reasonably well with the Onsager calculation. In the 
following we will assume for simplicity that ly = 1. The above fit is the only instance 
in which we have included A^ = 10 and A” = 15 and thus we believe that the obtained 
value of —B^ is an underestimation of the true critical point. 

• Magnetization and The Critical Exponent /3: The magnetization and the zero 
power are defined by 

m =< |TrM| > , y =< |TrMp > — < |rrM| . (3.6) 

Po =< {^TrMf > . (3.7) 
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Measurements of the magnetization m/N were performed near the extrapolated 
critical point B = —3.07 for C = 1.0 but inside the uniform ordered phase. These 
are then used to compute the critical exponent (3 by searching for a power law 
behavior. 

More precisely, we measure ln(m/A^) versus InA^ for each value of B very near and 
around B = —3.10, fit to a straight line in the range 20 < < 60 and compute the 

slope /3, then search for the flattest line, i.e. the smallest slope (3. This value marks 
the transition from the Ising phase to the disordered phase. Deeep inside the Ising 
phase the slope should approach the mean field value —1/4 which can be shown from 
the scaling behavior of the dominant configuration. After determining the critical 
value we then consider the value of B nearest to it but within the Ising phase and 
take the slope there to be the value of the critical exponent (3. In our example here, 
the flattest line occurs at B = —3.13 with slope —0.088(10) after which the slope 
becomes —0.109(11) at B = —3.14. The slope goes fast to the mean field value —0.25 
as we keep decreasing B. See figure (jH). Our measured value of the critical point i?* 
from the magnetization and of the critical exponent /3 are therefore 

5* = -3.13. (3.8) 

772 

In — = -0.109(11). In iV - 1.423(43) ^ P = -0.109(11). (3.9) 

• Susceptibility and Zero Power and The Critical Exponent 7: The measure¬ 
ment of the critical exponent 7 is quite delicate and will be done indirectly as follows. 
We rewrite the susceptibility in terms of the zero power and magnetization as 


X = < |TrM|2 > - < |rrM| 7^ 

= N^Po-m'^. (3.10) 

The critical exponent 7 in terms of the critical exponent j of Pq is then given by 

7 = 2 + 7 '. (3.11) 

By using the results shown on table ([I]) at B = —3.14, plotted on figure ([5]), we 
obtain the following exponents 

InPo = -0.352(10). In A^- 2.289(36) ^ 7 ' = -0.352(10). (3.12) 

Or equivalently 

InA^^p^ ^ 1.648(10). In AT - 2.289(36) ^ 7 = 1.648(10). (3.13) 

For consistency we can check that the second term in the susceptibility behaves using 
the result (13.9p as 

Inm^ = 1.782(22). In A^- 2.846(86) ^7 = 1.782(22). (3.14) 
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Our two measurements of the critical exponent 7 agree reasonably well with the 
Onsager values. 

If we try to fit the values of the susceptibility at its maximum shown in third column 
of table O, i.e. at the peak which keeps slowly moving with B, then we will obtain 
a very bad underestimate of the critical exponent 7 given by 

Inxmax = 0.515(08). In N - 0.652(30) ^ 7 = 0.515(08). (3.15) 

This in our mind is due in part to the dependence of Be on N and in another part is 
an indication of the critical slowing down problem showing up in the measurement 
of this second moment, i.e. the size of the fluctuations is observed to grow with 
N at the critical point but not at the correct rate indicated by the independent 
measurements of the zero moment and the magnetization. See figure ([5]). 

• Specific Heat and The Critical Exponent a: The sepcific heat is defined by 

C^=< > - < S >^ . (3.16) 

The critical point i?* as measured from the specific heat is identified by the inter¬ 
section point of the various curves with different N shown on figure ([3]). We get 

5* = -3.08. (3.17) 

This measurement is contrasted very favorably with the independent measurement 
obtained from the extrapolated value of Be shown in equation ()3.5p but should 
also be contrasted with the measurement obtained from the magnetization shown in 
equation (j3.8p . 

By using the results shown on table ([T]) at the critical point B = —3.08, plotted on 
hgure (l5|), we obtain the following exponent 
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In ^ = 0.024(9). In N - 0.623(31) ^ a = 0.024(9). 


(3.18) 


N 

B^,B^ = -3.07 

Ac 

(Cy)„B, = -3.08 

B <B^ = -3.13 


ioW<* 

10 

-3.20 

1.704(2) 

56.467(94) 

-3.14 

2.1776(12) 

6.256(6) 

15 

-3.16 

2.089(2) 

129.111(217) 

-3.14 

2.7750(14) 

4.315(4) 

20 

-3.14 

2.436(3) 

229.861(389) 

-3.14 

3.4423(15) 

3.571(2) 

25 

-3.13 

2.716(3) 

365.183(621) 

-3.14 

4.1759(16) 

3.220(2) 

30 

-3.12 

3.017(4) 

524.253(891) 

-3.14 

4.9772(16) 

3.042(2) 

36 

-3.11 

3.283(4) 

749.099(1267) 

-3.14 

5.8878(15) 

2.860(1) 

40 

-3.11 

3.515(4) 

941.139(1607) 

-3.14 

6.5134(14) 

2.782(1) 

50 

-3.10 

3.864(4) 

1461.597(2479) 

-3.14 

7.9250(12) 

2.576(1) 

60 

-3.10 

4.301(5) 

2144.929(3658) 

-3.14 

9.2021(11) 

2.388(1) 


Table 1: Measurements of the magnetization (m/iV)<*, the susceptibility x<*, via the zero 
power (Po)<*, and the specific heat (C^/iV^)* used to compute the critical exponents /9, 7 and 
a respectively. Here C = 1.0, the extrapolated critical point is B = —3.07, the critical point as 
intersection point of curves of specihc heat is B = —3.08, and the critical point as the flattest 
line of decrease of magnetization is B = —3.13. 

3.3 Free Propagator from Wigner Semicircle Law 

We can also measure the emergent geometry by measuring the free propagator of the 
theory. This will give us information on both the dimension and the metric since the free 
propagator is the inverse of the Laplacian A which fully encodes the underlying geometry 
in the sense of UM- This goes as follows [32]. 

A noncommutative phi-four on a d—dimensional noncommutative Euclidean spacetime 
Rg reads in position representation 

S = [ d^x{-di^di^ + ( 3 . 19 ) 

J 2 2 4 

The first step is to regularize this theory in terms of a finite Af—dimensional matrix d* and 
rewrite the theory in matrix representation. Then we diagonalize the matrix The mea¬ 
sure becomes f Hi d^iA^(^>) f dU where <l>i are the eigenvalues, A^(<1>) = 
is the Vandermonde determinant and dU is the Haar measure. The effective probability 
distribution of the eigenvalues <l>i can be determined uniquely from the behavior of the 
expectation values < f >. These objects clearly depend only on the eigen¬ 

values and are computed using a sharp UV cutoff A. If we are only interested in the 
eigenvalues of the scalar matrix <I> then the free theory A = 0 can be replaced by the 
effective matrix model |32] 

S = ^Trd>2. (3.20) 

«o 

This result can be traced to the fact that planar diagrams dominates over the non-planar 

ones in the limit A —?• 00 . This means in particular that the eigenvalues are distributed 
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according to the famous Wigner semi-circle law with being the largest eigenvalue, viz 


Pi't) = ~ > ““0 <t < -Kao- 

vron '' 


In the most important cases of d = 2 and d = 4: dimensions we have explicitly 
Q!o(m, A) = (A^ - n? ln(l ^)) , d = 4. 




(3.21) 


(3.22) 


1 

ao(m,, A) = — ln(l -|-^) , d = 2. (3.23) 

vr 

Obviously, dimension four is eliminated by the results of the critical exponents. In two 
dimensions the regulator A originates in only one of two possible noncommutative spaces 

m- 


1. Fuzzy Torus: As it turns the results on the fuzzy torus are different from those 
obtained using a sharp momentum cutoff due to the different behavior of the propa¬ 
gator for large momenta and as a consequence the resulting formula for Oq is different 
from the above equation (|3.23l) . We obtain instead 


ag(m,A) = 4 / 

Jo 


dh 


0 (27r)2 ^.(I - cosri)m2/2/2 ’ 


d = 2. 


(3.24) 


I here is the lattice spacing, the noncommutativity is quantized as 9 = NP/tt and 
the cutoff is 


A = 




(3.25) 


The above behavior can be easily excluded in our Monte Carlo data and by hindsight 
we know that this should be indeed so because the original multitrace approximation 
is relevant to the fuzzy sphere. 

2. Fuzzy Sphere: The fuzzy sphere = CPjy is the simplest of fuzzy projective 
spaces CPjy. In this case J\f = N + 1 and the scalar field <I> becomes an A^ x 
matrix (p given hy (p = y'^vr/iVad). In this case the cutoff is given in terms of the 
matrix size N and the radius R of the sphere by 


A = 


N 

'r' 


Also, in this case the mass parameters B and are related by 


m 


2 


b 


(3.26) 


(3.27) 


By using B = B/N^P and choosing a = Itt/N, so that ^ = cp, we obtain 


A^ _ 27r 

^2 ^/NB' 
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(3.28) 











We get then 


1 27r 

al{m, A) = - ln(l + -y^)- (3-29) 

In the limit B —> oo we get the one-cut = 2N/B of the Gaussian matrix model 
BTrlVP, viz B = 2M/a^. This can also be obtained by taking the limit B —> oo 
of the one-cut (deformed Wigner semicircle law) solution 

p{X) = -^i2CX^ + B + , 5^ = :^i-B + ^B^ + 12ArG) 

iVvr oC 

(3.30) 


of the quadratic matrix model BTrM‘^ + CTrM^. 

This result was also generalized in m- The eigenvalues distribution of a free scalar 
field theory on the fuzzy sphere with an arbitrary kinetic term, viz S = Tr{MlCM + 
BM‘^)/2, where /C(0) =0 and JC is diagonal in the basis of polarization tensors 
is always given by a Wigner semicircle law with a radius 


= <52 



N 


N-l 


/(B) = E 

1=0 


21 + 1 

K{l) + B' 


(3.31) 


Some Monte Carlo results are shown on figures Q and ©• These are obtained in 
Monte Carlo runs with 2^*^ thermalization steps and 2^® thermalized configurations where 
each two configurations are separated by 2^ Monte Carlo steps in order to reduce auto¬ 
correlation effects. We consider = 20 — 40, C = 0.05 — 0.35 and 5 = 0 — 5. 

It is not difficult to convince ourselves that the mass parameter B is precisely the 
mass squared in this regime. For each value of {N,C,B) we compute the eigenvalues 
distribution p{X) and fit it to the Wigner semicircle law (13.211) (see figure ([6])). We obtain 
thus a measurement of the radius of the Wigner semicircle law 6^ = cxq = We have 
checked carefully that in this regime the Wigner semicircle law is the appropriate behavior 
rather than the one-cut solution (|3.30p as evidenced by the first graph in figure (l6|). The 
measurement of the radii for various values of B is then plotted and compared with the 
expected theoretical behaviors (13.291) as well as with the B —00 behavior = 2N/B 
(see figure ([7|)). The agreement with (j3.29p is very reasonable with some deviation for 
small values of B as we approach the non-perturbative region where the uniform ordered 
phase appears at some 5 < 0. This discrepancy for small values of B is already seen on 
figure Q when we fit the distributions to the Wigner semicircle law. However, this effect 
is reduced as we decrease the value of C. 

In summary we conclude that we are indeed dealing with the geometry of the fuzzy 
sphere and, given hindsight, we know that this should be true. 


4 Generalization and Conclusion 

The emergence of geometry in the very early universe is a problem of fundamental 
importance to our understanding of quantum gravity and cosmology. In this letter, we 
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have proposed a novel scenario for the emergence of geometry in random multitrace matrix 
models which depend on a single hermitian matrix M with full unitary U{N) invariance 
and without any kinetic term. Thus, the model under consideration has no geometry 
a priori precisely because of the absence of a kinetic term. On the other hand, previous 
proposals of emergent geometry required the input of several matrices with some rotational 
symmetry group besides the U{N) gauge symmetry |35j . 

Our proposal consists in checking whether or not the uniform ordered phase is sus¬ 
tained by the multitrace matrix model under consideration. If yes, then the dimension 
of the underlying geometry, in the region of the phase diagram where the uniform or¬ 
dered phase is stable, can be inferred from the values of the critical exponents of the 
Ising phase transition. Whereas, the metric/Laplacian of this geometry can be inferred 
from the behavior of the free propagator encoded in the Wigner semicircle law behavior 
of the eigenvalues distribution of the matrix M in the weakly coupled regime. An explicit 
example is given in which the geometry of the fuzzy sphere emerges, with all the correct 
properties, in the phase diagram of a particular multitrace matrix model containing mul- 
titrace terms depending on the moments mi = TrM, m 2 = TrM‘^ and m 3 = TrM^ in a 
particular way m- 

This idea can be generalized in a straightforward way to all higher fuzzy projective 
spaces CP"" and fuzzy tori T"" by tuning appropriately the coefficients of the multitrace 
matrix model and/or including higher moments in the multitrace matrix model . 
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Figure 1: The phase diagram of noncommutative phi-four theory on the fuzzy sphere. In the 
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